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Abstract. We study the electronic band-structure and transport properties of a Cu02-plane within the 
three-band Hubbard model. The Dynamical Mean-Field Theory (DMFT) is used to solve the many particle 
problem. The calculations show that the optical gap A op t is given by excitations from the lower Hubbard 
band into the so called Zhang-Rice singlet band. The optical gap A op t turns out to be considerably smaller 
than the charge transfer energy A (A = e p — ea) for a typical set of parameters, which is in agreement 
with experiment. For the two-dimensional Cu02-plane we investigated the dependency of the shape of 
the Fermi surface on the different hopping parameters tcuO and too- A value foo/fcuO > leads to a 
Fermi surface surrounding the M point. An additional different static shift of the oxygen energies is also 
considered to calculate the electronic response due to a displacement of the oxygen atoms given by a 
frozen phonon. The density-density correlation for the oxygen orbitals is linear in doping for both hole 
and electron doping but shows a different temperature dependency in the two regimes. In the first case 
it is temperature independent and increases upon doping, which leads to an increasing electron-phonon 
coupling for the Bi 9 -mode in high-T c superconductors. 

PACS. 71.27.+a Strongly correlated electron systems - 71.30.+h Metal-insulator transitions and other 
electronic transitions - 74.25.Fy Transport properties (electric conductivity) - 74.25.Jb Electronic struc- 
ture - 74.25.Kc Phonons 



1 Introduction 

Since the discovery of high temperature superconductiv- 
ity by Bednorz and Miiller [0 it is a challenging prob- 
lem to find a theoretical description of this phenomenon. 
However, even in the normal state the electronic structure 
of the new compounds, which show superconductivity, is 
extremly hard to describe due to electronic correlations. 
A characteristic feature of the high-T c compounds is a 
unit cell with one or more Cu02-planes, which are re- 

| sponsible for the superconductivity §]. The most simple 
description of such a Cu02-plane is achieved within the 
one-band Hubbard model || , which considers one effective 

, 3d orbital of Cu in a tight-binding model in the presence 
of a local Coulomb repulsion Ud- This model describes 
well the properties of an insulator, because it reproduces 
the so called Mott-Hubbard metal-insulator transition Q . 
However, this simple model leads to a wrong description 
of the planes concerning the doping dependency of vari- 
ous properties, such as the asymmetric magnetic doping- 
temperature phase diagram |5j. 

The oxygen atoms in the Cu02-plane, respectively the 
2p-orbitals of these atoms introduce a further degree of 
freedom and the strongest hybridization takes place be- 
tween the Cu-3d x 2_ y 2 orbital and the 0-2p x / y orbitals Q. 
Therefore they form the lowest lying bonding state and the 



Send offprint requests to: 



highest lying anti-bonding state. Choosing the most sim- 
ple description of the Cu02-plane one is led to the hole 
picture and one has to consider the latter state, which is 
occupied with one hole. In order to describe the additional 
degrees of freedom, one has to solve at least a so called 
three-band Hubbard model or Emery model (?]] , with two 
additional oxygen orbitals 2p x and 2p y , which are included 
in a tight-binding manner. For the sake of simplicity all 
Coulomb energies but Ud are neglected. Within this model 
one is able to describe a doped charge transfer insula- 
tor, which gives a more realistic picture of the planes in 
a high-T c compound. The three-band model has recently 
been studied by Schmalian et al. using a generalized 
dynamical mean-field (DMFT) approach and our contri- 
bution is based on their approach. Watanabe et al. |J dis- 
cussed the metal-insulator transition in a two-band Hub- 
bard model in infinite dimensions and discovered the coex- 
istence of metallic and antiferromagnetic phases. This was 
also found by Maier et al. [jlO| , who additionally reported 
about an asymmetric magnetic doping-temperature phase 
diagram within our approach, which shows, that the anti- 
ferromagnetic phase is more stable upon electron doping 
than upon hole doping. 

In this paper the DMFT approach of Schmalian et al. 
H for the three-band Hubbard model is used. Several ex- 
tra features in two spatial dimensions were added, such 
as an oxygen-oxygen hopping process and a splitting of 
the oxygen energy levels. We discuss results for the one- 
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particle spectra of the orbitals and for the in-plane con- 
ductivity for hole and electron doped systems considering 
a semi-elliptic density of states. We further discuss, for 
the case of a two dimensional system, how an additional 
oxygen-oxygen hopping process affects the shape of the 
Fermi surface. Finally a static shift of the oxygen ener- 
gies due to a displacement given by a frozen phonon is 
considered in order to calculate the static density-density 
correlation, which describes the electronic response to this 
kind of distortion. A summary will conclude the paper. 

2 The three-band Hubbard model 

The starting point of our approach is a simplified Emery 
Hamiltonian JtJ describing the dynamics of holes in a doped 
Cu02-plane. The nearest neighbour hopping processes be- 
tween the 3d-orbital of the Cu-atom and the 2p x -/ 2p y - 
orbital of the O-atom (Uj v and tf iK j V ) are taken into ac- 
count. Sd, £ Px and e Py represent the energy levels of each 
orbital and fj, the chemical potential, in addition a local 
Coulomb energy Ud for the case of a doubly occupied d- 
orbital, which is responsible for the correlations. 

+ y ] (^ijV d i,aPjv,(T + h.C.) + y ] t iK j w p\ K a Pj v , a 

i 

(1) 

As mentioned earlier, this Hamiltonian is able to describe 
a doped charge transfer insulator, which is characterized 
by the charge transfer energy A = e p — Ed- The relation 
U > A ensures the system of being in the charge trans- 
fer regime pi. This relation typically holds for parame- 
ters obtained from a first principles calculation (see table 
[l]). The gauge invariance of the Hamiltonian allows us to 

Table 1. Parameters for a three-band model (in eV) calculated 
with a constrained first principles calculation for La2Cu04 
done by Hybertsen et al. 113]. 



A t t' 


u d 


U P 


U p d 


Upp 


3.6 1.3 0.65 


10.5 


4 


1.2 






choose the phase of the wave functions freely. We use the 
phase convention shown in figure [l], which determines all 
phases of the considered hopping processes. The transfor- 
mation into k-space leads us to the following Hamiltonian 
(for simplicity we show the non-interacting part only): 

h = J2 (4,.pL,.pU) KH > ( 2 ) 

k,<T \Pyk<? J 




Fig. 1. Used phase convention, which determines the phases 
of all hopping processes. 

with the matrix /i(k): 

(Ed — M 2i£sin(kr x ) 2itsin(kr y ) \ 

— 2iisin(kr x ) e Px — /J, — it' sin(kr x ) sin(kr y ) I 

— 2itsin(kr y ) — it' sin(kr x ) sin(kr y ) e Vy — fj, J 

(3) 

This uncorrelated Hamiltonian can be diagonalized easily. 
In order to treat the strong Coulomb repulsion Ud-, we use 
the dynamical mean-field theory, where the lattice prob- 
lem is mapped onto an effective Anderson impurity model 
[ p"3| . This model describes an impurity and its hybridiza- 
tions with conduction electrons (holes in our case). The 
map ping onto this model is exact in infinite dimensions 
[y_5| and means the restriction to local selfenergy diagrams 
in finite dimensions. 

As first step of our derivation of the cluster/impurity we 
apply a unitary transformation to the Hamiltonian (2): 

(e A - fi -2t7 k \ / d k<7 \ 

k,<X \ t'y. £j5k _ M/ \Pk<T J 

(4) 

with 

£ P k = -^{sin(kr x ) 2 £ Px +sin(kr y ) 2 £p H - 
7k 

- 8t' sin(kr x ) 2 sin(kr y ) 2 }, (5) 
£pk = -^{sm(kr y ) 2 £ Px + sin(kr x ) 2 e Py + 

7k 

+ 8i'sin(kr x ) 2 sin(kr y ) 2 }, (6) 

, sin(kr x ) sin(kr v ) 

j = — v x % v yJ x 

{(£ P , " £ Py ) + 4i'(sin(kr x ) 2 - sin(kr y ) 2 )} (7) 

and 

7 2 =sin(kr x ) 2 +sin(kr y ) 2 . (8) 
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This transformation leads to new orbitals p and p. d hy- 
bridizes with p but not with p. The case of s Pal = e Py = e p 
and t' = leads to a dispersionless p-band at e p — /i and 
the p-orbital is decoupled from the rest of the system. 
This procedure leads to a well-defined description of local 
orbitals at different lattice cells, which was suggested by 
Valenti and Gros O]. Within this formulation the prob- 
lem how to distribute the p-orbitals, which orginally are 
located in between different Cw-sites, to a particular unit 
cell was solved. In order to construct a sensible DMFT we 
here follow a path, which identifies the effective 'impurity' 
needed to set up the DMFT equations directly from the 
d-part of the Green's function. Due to the form (4) of the 
Hamiltonian the general structure for G^(z) within the 
DMFT is: 



z — £ P k + A* 



t' 



(9) 



The sum over k will in general lead to a complicated struc- 
ture in the denominator. However, since the sclfencrgy 
Ed(z) is k-independent we may cast it into the form: 



of the effective impurity is defined by: 

a a 
a 

+ y^(^k Ptc^a + h.C.) + ^ £k C k ,<j c k,<r , 
k,tT k.er 



(12) 



with a hybridization Vk of the p-orbitals to the effec- 
tive medium described by the dispersion e^. The effective 
medium enters the solution of the impurity problem only 
by the hybridization function 



4M = £E 



k 



(13) 



of the p-states, which has to be determined selfconsistently 
from equation (10). 

Thus equation (10) is the central point of our selfcon- 
sistency scheme. After calculating the local Green's func- 
tion by the sum of k one can extract the hybridization 
function A(z), which is used to solve the impurity prob- 
lem with help of an extended NC A scheme g . From the 
impurity Green's function for the (i-orbital we extract the 
selfenergy 2J(z), which is then used in the calculation of 
the lattice Green's function in equation (10). 



k 

z - e d + n - E d (z) ■ -—^ , 

z - £ P oo + M - A(z) 

(10) 

where t 00 and £ p oo are Fourier-transformed terms of — 2£7k 
and £ p k- The effective medium's function A(z) is defined 
through equation (10) and incorporates all bandstructure 
effects due to the dispersion and hybridizations. Note that 
this definition is far from being unique! In fact we may 
choose for example 



3 The effective two-band model 

Let us begin the discussion of the properties of the three- 
band Hubbard model by considering the standard case 
£ p = e Plc = e Py and t' — 0. Then the non- hybridizing 
p-orbital decouples from the rest of the system leading 
to a two-band problem. This problem can in principle be 
generalized to D dimensions, but it is impossible to find 
a uniform scaling to obtain a nontrivial limit D — > oo 
|L4| . The DMFT used here therefore always leads to an 
approximate treatment of the model. In order to simplify 
numerical calculations the semi-elliptic density of states is 
used, i.e.: 



1 J2 Gi d a (z) l{z-e d + v- S d (z) - A(z)}- 1 (11) Gd(z) = 1 £ ^ 



k 



z - (e d - fi) - E{z) - 



z-(e p - fi) 

Lt 2 (l - e) ' 
z- (e p - /x) 



(14) 



f j p (e) 

as another possibility. In this case the equation defining — / de ^ 2 

the DMFT would be formally the same as for the one- z — (ed — ji) — E{z) \ 

band Hubbard model. However, due to the singular struc- 
ture induced by a d-p-hopping equation (10) turns out to 
be numerically more convenient: For A(z) = the level 

structure of the d-p complex is already included from the with Po(e) = ^ V 1 - £ 2 . The Green's function for the p- 
outset and A{z) is rather smooth. On the other hand the orbital is obtained, if one exchanges the energies e d + E{z) 
at the first glance more natural choice (11) would result in anc ^ £ v- 
a A(z) with a highly singular behaviour reflecting the ex- 
istence of unrenormalized p-states. Thus the Hamiltonian 
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3.1 One-particle spectra 

Figures and ||b show the one-particle spectra for the d- 
orbital and the hybridizing p-orbital for an electron and a 
hole doped system, respectively, at room temperature, i.e. 
at a reciprocal temperature of (3 = 40 eV~ x . They show 
in terms of the hole picture from the left to the right: the 
lower Hubbard band, the Zhang-Rice band mostly due 
to exitations into a two particle singlet state, the p-band 
and the upper Hubbard band, which has a double peak 
structure due to the complex excitation capabilities of the 
cluster under consideration. In both figures typical particle 
hole excitations are marked by arrows, w hich play a role in 
the dynamical conductivity (see section 3.2). In addition 



to those expected structures a sharp resonance appears 
near the Fermi energy in both doping regimes. From the 
DMFT of the one-band Hubbard model it is well known 
JT(| , that this resonance is connected to a local Kondo-like 
screening. It thus is frequently termed Abrikosov-Suhl res- 
onance to stress the similarity of its physics to the Kondo 
effect. For a detailed description of the doping and tem- 
perature dependencies of the spectral weigths see ref. M . 



3.2 The in-plane conductivity 

In order to calculate the in-plane conductivity in linear 
response theory using Kubo's formula, one has to derive 
the current-current correlation function. The current den- 
sity operator for the three-band model ( see equation (1)) 
is given by the time derivative of the polarization opera- 
tor using Heisenberg's equation of motion. This procedure 
leads to the expression: 



2et 



EE: 

v k.i 



v v sin(kr„)(dJ tj(T p J ,k, ( T + P v ^ a dk, a ), 



with the vector r„ connecting d- and p-sites. 

In the following the real part of the in-plane conduc- 
tivity in x-direction is examined. The current-current cor- 
relation function is calculated within the orbital represen- 
tation of equation (4). Since the particle-hole excitations 
with the non- hybridizing p-states are playing a role only at 
high energies, contributions with p-orbitals are neglected. 
Thus we obtain: 



1 .aet 2 1 ^ 



sin(kr a; ) 4 + sin(kr a ) 4 



/(e) -f(e + co) 
de i7(e,w,k) 



with 



iT(e,w,k) 



E p° 



w,k)p (e,k) 




7.5 10.0 12.5 

0) in [eV] 




Fig. 2. One-patricle spectra for the d-orbital and the hybridiz- 
ing p-orbital in a electron doped system (a) and in a hole doped 
system (b) for energies in units of the hoppping parameter t. 
The parameter set was f3 = 40.0 eV' 1 , A = 3.6 eV, t = 1.0 eV 
and U d = 7.2 eV. 



and spectral function p a of Green's functions with the 
following orbital combinations: 

(a, /3) e { (ppt , drft ) t (dd!,pp*), (pit , pS ) , (dpt , d p1 )} . 



Due to the locality of the vertex kernel and the antisym- 
metric current operator vertex corrections vanish analo- 
gously to calculations for the one-band Hubbard model 
[y_7|. Finally the k-summation is carried out in an aver- 
aged way by integration over the pure density of states: 



OO rOC 



dede p D (e) U(e,(j,e). 



OO J —OO 



Typical particle-hole excitations for electron and hole 
doped systems are denoted by the vertical arrows in fig- 
ure ^|a and ||b. These different exictations can in principle 
all be resolved in the graph of the conductivity. However, 
for sake of simplicity only the low energy scale is shown 
in figures and [3^3. The in-plane conductivity was cal- 
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Fig. 3. Real part of the conductivity for several electron doped 
(a) and several hole doped systems (b) at f3 = 40 eV -1 . The 
inset shows the temperature dependency for a system with the 
doping level x = — 4 %. 



4 Results for a quadratic Cu0 2 -system 

In this section we present results for a two-dimensional 
tight-binding model. At a first glance the use of the DMFT 
as an approximation may seem to be particularly crude. 
However, as long as the system is not too close to a phase 
transition the choice of a local selfenergy appears to cap- 
ture the most important dynamics, which is seen if one 
compares the results from a local approximation with those 
of a Monte Carlo calculation in the two dimensional one- 
band Hubbard model p^j . 

4.1 Bandstructure and Fermi surface 

Before we discuss the bandstructure for the correlated 
electron system, let us look at the uncorrelated tight- 
binding bandstructure, which is shown in figure §. In the 
following the name of the bands is given by their main 
orbital character. Here the band with the lowest energy 
(hole picture) has primarily d-character. The upper bands 
have p-character. For hole-doping the Fermi energy lies in 
the upper bands, for electron-doping in the lower band. 
The d-band has a minimum at the M-point and a sad- 
dle point at the X-point. With increasing value of t'/t 
the distance between the minimum and the saddle point 
increases. This means that the region of electron doping, 
where one can find a Fermi surface surrounding the M- 
point by adjusting the chemical potential /i, is also increas- 
ing. In the hole doped case, however, a positive value of 
t'/t leads to a Fermi surface centered around r produced 
by the hybridizing p-band and one around M by the non- 
hybridizing p-band, since the non-hybridizing p-band is 
pushed towards the hybridizing p-band with an increas- 
ing value of t'/t. On the other hand, a value t'/t < 

6 



culated for several doping parameters x for a reciprocal 
temperature (3 — 40 eV~ x . The temperature dependency 
for constant doping x is shown in the inset of the figure ||a. 
Both electron and hole doped systems show an increasing 
conductivity upon doping in the mid-infrared region at 
0.5 eV. The far-infrared part of the conductivity shows a 
pronounced increase due to the growth of spectral weight 
of the Abrikosov-Suhl resonance upon doping. In the in- 
set of figure ||a the strong temperature dependency of the 
low energy excitations is shown. The broad peak at 3 eV 
shown in the inset corresponds to transitions between the 
lower Hubbard band and the Zhang-Rice band. The en- 
ergy separation of these two bands, which describes the 
optic energy gap A opt , is smaller than the bare charge 
transfer energy A = e p — Ed- For example, for the present 
set of parameters we find A opt rs 2 eV, which is in good 
agreement with meassurements by Uchida et al. [18l for 
La2- x Sr x Cu04 and Nd2- x Ce x Cu04_ y . 




r m x r 

Fig. 4. Tight-binding bandstructure for several values of t'/t. 
e d = eV, e p = 3.6 eV, t = 1 eV and t' = 0, ±0.3 eV. 

would push down the non-hybridizing p-band, leading to 
a dispersion minimum at the M -point, as it is observed 
experimentally pof . Thus the requirement of a Fermi sur- 
face surrounding the M-point for both electron and hole 
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doping means that one either has to assume different signs 
for t' in the two cases of doping or look for a more subtle 
mechanism leading to the observed physics for a fixed sign 
oft'. 

In the following we want to argue that in order to solve 
the puzzle of the correct choice of the parameter t' , one 
has to consider electronic correlations. In order to make 
contact with the free bandstructure we plot the total spec- 
tral weight A(k, u) in a density-plot in the cj-k-plane. The 
dark regions refer to high spectral weight. In figure ^|a and 
||b the low energy part of the resulting k-dependent total 
spectral weight for the electron and hole doped CuCV 
plane is shown. These plots show the lower Hubbard band, 
the Abrikosov-Suhl resonance and the Zhang-Rice band. 
These bands exhibit the typical dispersion of the tight- 
binding d-band. On the other hand the p-bands and the 




r m x r 




-a A(k,co) k 

r M r 

Fig. 5. Total spectral weight in a density-plot in the cj-k-plane 
for special k-points for f3 = 40 eF~\ t = 1 eV, t' = 0.3 eV, 
A = 3.6 eV, U d = 7.2 eV and x = -8.5% (a),i = +8.5% (b). 

upper Hubbard band at higher energies (not shown in the 
figure) have the dispersion of their uncorrelated pendants. 
For both types of doping there occurs a quasi-particle 
Abrikosov-Suhl resonance at the Fermi energy. It is impor- 
tant to note, that the dispersion of this band (which has 
d-character) completly determines the shape of the Fermi 



surface. In figure r| we show, that for constant doping pa- 



M 




r x 

(a) t' = -0.2t 

M 




r x 

(b) t' = +0.2t 

Fig. 6. Total spectral weight at the Fermi energy in a density- 
plot for two values oft'/t in the irreducible part of the Brillouin 
zone. The doping is x = +5% in both cases, the remaining 
parameters are the same as in figure ^. One can see that the 
choice of t'/t determines the shape of the Fermi surface. 

rameter x — +5%, one can change the shape of the Fermi 
surface by varying the value of t'/t, since it directly influ- 
ences the dispersion of the band. From the experimental 
fact that the Fermi surface encloses the M-point we are 
led to a value t'/t > 0, which is in agreement with those 
anticipated for high-T c compounds such as Nd2- x CeCu04 
@ or Bi 2 Sr 2 CaCu 2 8+x f§. 

4.2 Electronic Response to a frozen phonon excitation 

In this section we investigate the response of the oxygen 
occupation on a level splitting of the p x - and p y -states. 
Such a level splitting may be realized by a different dis- 
placement in z-direction of the x- and y-oxygen in the 
Cu02-plane. For simplicity we treat a symmetric energy- 
level splitting, i.e. one has to modify equation (7) to (9) 
with e v i = e v ± Se v . With this modification we can look 
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at the static density-density correlation function in linear 
response approximation: 

8{n x ) - 8 (n v ) 



{{8fl X y,5h X y))^ = Q) 



(15) 



This quantity can be used to examine the electron-phonon 
coupling for the -Bi ff -phonon in high-T c -compounds. Of 
course it is not possible to study the dynamical aspects of 
the coupling, but it is a rough meassure for the renormal- 
ization of the phonon-frequency due to this coupling. 

Up to Se p = 0.075 eV the correlation function was 
found to be independent on Se p for all considered sets of 
parameters. The static susceptibility was therfore inves- 
tigated for a constant value of the level splitting, 5e p — 
0.05 eV. In figure [?] we show the density-density correla- 
tion for t'/t — 0.5 for several temperatures as a function 
of the doping x. The dashed line marks the doping which 



O p= 10.0 1/eV 
□ 40.0 1/eV 
O 100.0 1/eV 



A 
A 



S 

to 

V 
V 




doping level x in [%] 

Fig. 7. Static density-density correlation at three different 
reciprocal temperatures /3 for Se v = 0.1 eV, t'/t = 0.5, 
U d = 7.2 eV and A = 3.6 eV. 



is found, when the chemical potential /i is in the gap. Note 
that the crossing point is not exactly at x = 0. Thus the 
doping x cannot be identified with the doping level of a 
real Cu02-plane. This shortcoming originates from a nu- 
merical failure of the NCA and is enhanced by a large 
value of t'/t — 0.5. This effect can be neglected in the 
previous section where for t' = the model doping 
level x for the insulating case was less than 1%. For clar- 
ity we term the region left to the dashed line in figure [7] 
as electron doping and the one right as hole doping. In 
both regimes the density response is a linear function of 
the doping x. This holds for all investigated sets of pa- 
rameters. In the hole doped regime the density-density 
correlation shows in addition no effect upon changing the 
temperature and it can be described by a straight line with 
a slope d((Sn xy , 8h xy )) ( u =o) /dx « —6 eV^ 1 . For electron 
doped systems the slope increases with decreasing temper- 
ature. In figure |8| the hole-density response of p-orbitals is 
shown together with occupation of p- and d-states for dif- 
ferent doping. It can be seen that the desity response is 
strong in the hole doped regime (x > 3%), when also the 
number of holes in the p-states is large. 




0.0 3.0 
doping level x in [%] 

Fig. 8. Comparison of the static density-density correlation 
with the occupancies of the different orbitals for 5e p — 0.1 eV, 
t'/t = 0.5, U d = 7.2 eV, A = 3.6 eV and f3 = 40 eV -1 . 



We have also calculated the hole-density response func- 
tion for d-states at (71", 7r). This function is rather indepen- 
dent of the Coulomb parameter Ud, but depends on the 
occupancy of the (i-orbitals. For a calculation with t' = 
the strength of the d-response varies by a factor 3 to 4 
compared to the p-response due to the larger number of 
occupied d-orbitals. 



5 Summary 

In this paper we presented results for the electronic prop- 
erties of the three-band Hubbard model obtained with the 
help of the Dynamical Mean-Field Theory. For the calcu- 
lations we used two different model density of states: a 
simple elliptic density of states to calculate the local spec- 
tral density and the optical conductivity, a two dimen- 
sional tight-binding density to calculate the details of the 
bandstructure and a static response function. 

The spectral density consists of a lower Hubbard band 
with d-character, transitions into a Zhang-Rice band of 
p-d-singlet states, binding and non-binding p-bands, and 
a upper Hubbard band. For electron doping the Fermi 
level is in the lower Hubbard band, for hole doping in the 
Zhang-Rice band. The energy separation between these 
two bands, which is the optical gap A opt typically seen 
in conductivity meassurements, is smaller than the bare 
energy separation of p- and d-states, A — s p —Sd- This ten- 
dency as well as the calculated value of A opt qualitatively 
agree with experimental observed data. 

At low energies the optical conductivity shows a nar- 
row Drude-like peak on a broad background as experimen- 
tally observed. The narrow Drude peak is due to tran- 
sitions inside the Kondo-like feature close to the Fermi 
energy. 

We studied in detail the influence of an additional 
tigth-binding coupling t' between oxygen states on the 
bandstructure. It turned out that in the presence of cor- 
relations the variation of the Fermi surface is consider- 
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ably changed compared to the non-interacting system. 
The most important aspect is that one finds similar struc- 
tures of the Fermi surface for both electron and hole doped 
systems. Together with the experimental fact that the 
Fermi surface encloses the M-point for both kind of doped 
systems we conclude that t' jt must be positive. 

We also investigated the response of the occupation- 
difference of p-states on the energy splitting of the two 
different p-states in the unit cell. Such a response function 
is a meassure for the electron-phonon coupling of special 
phonons. Naturally this density response increases with 
the occupation number of p-holes, i.e. it is larger if the 
Fermi energy is in the Zhang-Rice band than in the lower 
Hubbard band. 
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